Microscale swimming: The molecular dynamics approach 
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The self-propelled motion of microscopic bodies immersed in a fluid medium is studied using 
molecular dynamics simulation. The advantage of the atomistic approach is that the detailed level 
of description allows complete freedom in specifying the swimmer design and its coupling with the 
surrounding fluid. A series of two-dimensional swimming bodies employing a variety of propulsion 
mechanisms - motivated by biological and microrobotic designs - is investigated, including the use 
of moving limbs, changing body shapes and fluid jets. The swimming efficiency and the nature of the 
induced, time-dependent flow fields are found to differ widely among body designs and propulsion 
mechanisms. 

PACS numbers: 47.63.Gd, 47.11. Mn, 87.19.St 



I. INTRODUCTION 

There is a growing demand - driven by nanotechnol- 
ogy - to understand and utilize nanoscale hydrodynam- 
ical phenomena, a consequence of which is the need to 
achieve a capability for detailed modeling of the underly- 
ing physical processes occurring at submicroscopic scales. 
There are numerous instances where nature itself pro- 
vides examples of successful design at this level, a regime 
where behavior runs contrary to the familiar macroscopic 
world; emulating some of the multitude of biological ar- 
chitectures and mechanisms may, therefore, provide in- 
sight into the requirements for optimized engineering de- 
sign. 

Microrobotics has long captured the imagination [1], 
with practical implementations only now becoming fea- 
sible. One category of microrobot with potential medical 
applications employs self-propelled swimming; however, 
long before such devices were envisaged, studies of the 
dynamics of microscopic life forms, and the theoretical 
analysis of hydrodynamics at appropriately low Reynolds 
numbers. Re ~ 10~^, made it clear that the prevailing 
fluid environments were unlike those encountered in more 
familiar circumstances (where Re ^ 1) but rather re- 
sembled slow swimming in highly viscous treacle. This 
is known as Stokes (or creeping) flow, and responsibil- 
ity for the counterintuitive behavior lies in the absence 
of inertia from the dynamics [2]. Microscopic creatures 
swim by employing mechanisms appropriate to these con- 
ditions, such as rotating flagella, beating cilia, body de- 
formation, and longitudinal or transverse surface waves 
0, 0, HI- While microscale robots may need to travel 
faster than their target biological organisms in the fluid 
medium, their motion relative to the fluid would still be 
at low Re. 

Low- Re flow is amenable to theoretical analysis be- 
cause the limiting form of the Navier-Stokes equation 
is linear and time independent [2]; in this limit, the 
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motion of self-propelled bodies can be studied using 
suitable approximations, perturbation methods and sim- 
ple models, e.g., |, 0, i, i, 0, [O, E, [H. Micro- 
scopic swimmers have also been synthesized experimen- 
tally [14]. Swimming has been studied numerically, based 
on a continuum fluid (with grid discretization) and the 
immersed-boundary method for dealing with the inter- 
face between the fluid and an elastic body [l5[, and 
mesoscopically described fluids utilizing, e.g. , the lattice- 
Boltzmann method, have also been used [16]. 

The present study of microswimmers is based on 
molecular dynamics (MD) simulation [13], in which the 
fluid itself consists of discrete atoms that interact directly 
with the swimmer. The MD approach is flexible in re- 
gard to the level of structural detail that can be incorpo- 
rated and, unlike other approaches, hydrodynamic cor- 
relations emerge naturally. A selection of self-propelled 
bodies that swim using a variety of mechanisms is in- 
vestigated; a particularly important aspect of the results 
is a comparison of the efficiency of the alternative body 
designs and propulsion techniques. 

The MD approach is motivated and justified by its suc- 
cess with other hydrodynamical problems, e.g., hexago- 
nal flow patterns in Rayleigh-Benard convection [18] and 
toroidal rolls in Taylor-Couette flow [H, [2^ , the latter 
even in quantitative agreement with continuum theory. 
Although large driving forces and gradients are required 
to compensate for small size in order to exceed the in- 
stability threshold, MD demonstrates continuum-like be- 
havior in fluid layers a mere 30 atomic diameters thick. 
Thus the scaling implicit in dynamic similarity extends 
down to the MD regime, although there is an eventual 
lower size limit [21]. 



II. METHODOLOGY 

The swimmer is immersed in a two-dimensional (2D) 
fluid (where computation requirements are reduced sub- 
stantially and flow visualization simplified compared to 
3D) confined to a rectangular container. The body and 
limbs of the swimmer are constructed from soft-disk 
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atoms. Each atom is placed at a particular location rela- 
tive to a reference coordinate frame attached to the body; 
the locations are either fixed in the frame if they form the 
body, or follow predetermined periodic paths if they be- 
long to a limb or other moving component of the propul- 
sion mechanism. The atoms interact via a short-range, 
repulsive interaction, U{r) = 4e[(cr/r)-'^^ — (a/r)^], with 
range r < Tc = 2^^^a. The atom spacing in the swim- 
mer is l.lcr, producing a rough boundary to oppose fluid 
slip; while a hindrance at higher Re where drag reduction 
is important, skin drag is necessary for certain kinds of 
low-Re propulsion. Swimmer mass is proportional to dis- 
placed fluid area (excluding limbs); while motion is not 
totally inertia free, it approaches this limit at low speeds. 
The container boundaries are made from flxed atoms, 
spaced to produce impenetrable, nonslip walls (periodic 
boundaries are unsuitable as they would allow bulk flow). 
Finally, there are the moving atoms of the fluid; fluid 
properties such as viscosity are intrinsic to the model 
and depend on the intermolecular forces, density, and 
temperature (in contrast to a continuum representation 
where viscosity is normally an assigned parameter). 

Interactions among fluid atoms, as well as between 
fluid atoms and the swimmer components and walls are 
computed at each MD timestep; the total force acting on 
the swimmer (and torque, if needed) is evaluated from 
contributions of all fluid atoms within range Tc of any of 
the swimmer atoms. Standard MD techniques [13] are 
used in the force computations (for larger systems the 
task can be parallelized), the equations of motion of the 
fluid atoms and swimmer center-of-mass (cm.) are in- 
tegrated using a leapfrog algorithm, and a thermostat is 
applied to prevent viscous heating. 

A broad range of swimmer designs employing differ- 
ent propulsion modes are considered; they are described 
below, and several are shown in Fig. [H Some designs 
aim to emulate biological life forms manifest at various 
size scales 0, [H, Hll, while others implement theoreti- 
cal models (organisms that appear biologically infeasible 
might be realizable artiflcially). Most designs are based 
on a cylinder - where the term 'cylinder' refers to the 
2D projection of a cigar (or spherocylinder) - with the 
same frontal area and body length to allow comparison of 
relative performance. Each swimmer executes closed pe- 
riodic cycles in the appropriate internal-coordinate space; 
most cycles are irreversible, an essential requirement for 
propulsion in the Stokes regime [3]. Swimmer names are 
indicative of shape or propulsion mechanism. 

Biflaps: two pairs of lateral flaps near the front and 
rear rotate from the forward direction through 90° and 
retract. Cigar: two curved paddles emerge laterally near 
the front, slide along the body to near the rear and 
retract - multiple paddles would resemble a transverse 
wave Flagellum: a flexible rear-pointing tail along 
which a transverse, two-cycle sinusoidal wave propagates 
- the 2D equivalent of a rotating helix. Flaps: a single 
pair of lateral flaps near the front. Jet: an open cavity at 
the rear into which paddles emerge from the inner walls 






FIG. 1: A selection of swimmer designs: biflaps, cigar, flagel- 
lum, jet, legs and snake; flow fields are shown. 

and slide rearwards acting as a piston to expel the fiuid, 
the paddles retract and there is a recovery period for re- 
filling the cavity - an approximation of squid propulsion. 
Legs: a pair of fiaps emerges laterally at the rear, rotates 
backwards through 90° and retracts. Slime: a circular 
body that extends a thin hollow rod (to a distance equal 
to the initial radius) at the end of which a new circu- 
lar body grows (at constant areal rate) while the initial 
body contracts (total area is fixed), and the rod finally 
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retracts into the new body - a related system is treated 
theoretically in [12]. 

Snake: based on a cylinder whose central axis is dis- 
placed laterally to follow a two-cycle traveling sine wave. 
Tail: a rigid, rear-pointing oar oscillates sinusoidally over 
a 36° range; this mechanism, unlike the others, is time- 
reversible, and under Stokes conditions swimming is im- 
possible [3] . Track: the atoms forming the sidewalls slide 
along the body, emerging from behind the front cap and 
disappearing into the rear, and propulsion is due to the 
skin drag of the rough boundary; this is an approxima- 
tion of a membrane extruded at the front and absorbed at 
the rear, or the longitudinal surface waves of cilia 0, Hf. 
Trilinear: three linked circular bodies, joined linearly by 
rods whose lengths vary with time in an irreversible cycle 
jllj . Tube: hollow cylinder open to the fluid at both ends 
with the same frontal area as the other cylindrical bodies; 
two paddles emerge inwards, slide to the rear expelling 
fluid, and retract - a 'peristaltic' mechanism. 

The swimmer cylinder length is L5 = 75; reduced MD 
units are used, where the length unit is cr = 0.34 nm, 
so that the actual length is 26 nm. Although this is 
several orders of magnitude less than even the smallest 
of organisms, the MD flow studies cited earlier indicate 
that this should not adversely affect the results. Other 
swimmer dimensions are proportional to L5; the radius 
is = I/5/6 (also the maximum radius of slime), so the 
total body length is 4L5/3 = 100, flap and leg lengths are 
L5/6, tail odii length L5/3, snake amplitude L5/IO; these 
and other design elements represent reasonable (though 
arbitrary) choices. Careful design is required to prevent 
fluid atoms becoming trapped in a corner by a moving 
limb; since only the time-dependent limb position is spec- 
ified, this would cause numerical instability. Swimmer 
components whose size varies (such as the growing flap) 
accomplish this by altering the atom overlap. The swim- 
mer cm. is constrained to move in a straight line and 
body rotation is not allowed, otherwise travel direction 
would be subject to random change (Brownian motion), 
complicating the analysis. 

A large container is required to help dissipate induced 
flows and ensure the swimmer remains far from the walls; 
container size is 835 x 338, and for fluid density 0.7 the 
total number of atoms is approximately 190 000. Runs 
begin with the swimmer at rest (except for towed bod- 
ies), positioned several body lengths from the rear wall; 
the fluid atoms are assigned random velocities with a rms 
value of unity, and arranged on a triangular lattice (out- 
side the swimmer). Run length is 610000 timesteps (the 
first 10 000 are excluded from analysis) of size 0.005 (MD 
time unit ^Jma^e = 2.2 ps, where m is the atom mass). 

The rate of change of the internal swimmer state is 
governed by the drive speed Vd] this corresponds to the 
(fixed or highest) speed of a sliding component or the tip 
of a rotating limb in the swimmer reference frame. The 
value of Vd determines the cycle period and, after allow- 
ing for interaction with the fluid, the average swim speed 
and power dissipation; Vd must be relatively small to min- 



imize density fluctuations in the (compressible) fluid and 
local shear rates, but sufficiently large that runs are of 
reasonable duration. 

Two externally towed bodies are included for compari- 
son, moving at fixed speed irrespective of the fluctuating 
fluid forces. Cylinder: the basic body design used for 
most swimmers. Disk: a circular body with radius 
but without the extra skin drag of the cylinder body. 



III. RESULTS 

Examples of the structured, time- varying flows that 
develop (for Vd = 0.6) appear in Fig. [T] (fluid atoms are 
not shown); swim direction is to the right. The spatially 
coarse-grained flow is evaluated over a rectangular grid 
attached to the swimmer, with velocities measured in a 
stationary frame; the swim cycle is divided into 16 phase 
segments and each is averaged separately over the run. 
Only a single segment is shown for each swimmer, but the 
phase chosen is where the flow best reflects the nature of 
the propulsion mechanism. Only the fluid region close to 
the swimmer is included to ensure adequate resolution; 
furthermore, the longest arrow for each swimmer corre- 
sponds to the highest flow speed anywhere in its cycle, 
and may not appear in the selected phase. 

Fig. [2] shows the average power (P) required to over- 
come fluid resistance and maintain a mean swim speed 
{vs) for several of the swimmers; the points on each per- 
formance curve are for different Vd or towing speed (in 
most cases 0.1, 0.2, 0.4 and 0.6, but reduced by factors 
of 2, 4 and 10 for slime, track and the towed bodies, re- 
spectively). P is the total work per unit time performed 
by swimmer atoms (both in the body and the propul- 
sion mechanism) against the forces exerted by fluid atoms 
within range Tc- {vs) is small relative to the fluid ther- 
mal velocity. The results are consistent with (vs) being 
proportional to Vd^ a consequence of the linearity of the 
flow problem. 

Higher efficiency implies reduced power for a given 
swim speed; the ratio {P)/{vs) provides a measure of 
performance - the lower the better. Efficiency differs 
considerably among swimmer designs, and for a given Vd^ 
some are faster (possibly utilizing more power) than oth- 
ers. Relative to the towed reference cylinder (which pro- 
vides a measure of the useful work done), track achieves 
a comparatively high efficiency of about 17%, as well as 
an effective coupling of the body wall to the fluid since 
(vs) ~ Vd/3, while for cigar efficiency is 2%, and only 
1% for tube; such results lie in the expected range [8(. 
The remaining swimmers either have lower efficiency or 
move even more slowly (briefly, for Vd = 0.6, flaps, flag- 
ellum and tail have ('Ug) = 0.04 ± 0.01 at (P) ^ 5, jet has 
(vs) = 0.04 at (P) = 11, while trilinear is even slower 
and less efficient). It should be noted that these compar- 
isons are applicable to the specific swimmers described 
here; the designs have not been optimized for maximum 
performance. 
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speed 

FIG. 2: Cycle- averaged power as a function of swim speed 
for different swimmers and drive speeds (MD units); separate 
graphs (with different scales) show (a) the most efficient of 
the swimmers and the towed bodies, and (b) swimmers of 
moderate efficiency. 

Establishing the relevance of the MD approach requires 
an estimate of Re. For the midrange speed {vs) = 0.04, 
and a measured kinematic viscosity v = 2.1 (from an MD 
analysis of 2D pipe flow - not shown), Re = 2Ri){vs) /i^ ^ 
0.5, a value where inertial effects can be disregarded [2^. 
Since the limb speeds are higher, inertial effects might be 
present at smaller length scales (explaining the motion 
of tail). The towed-body data can be fit to (P) = Cv^^ 
the functional form of simple Stokes drag, with Cdisk = 
11.5 and C cylinder = 18 (< 0.003 deviation); the latter 
is larger due to skin friction of the elongated body, in 
addition to pressure drag common to both. 

Fig.[3]shows the dependence of the instantaneous speed 
and power, Vs and P, on phase angle, for Vd = 0.4; 
the results are averaged over the corresponding phases 
of all cycles in the run (15-35, depending on swimmer). 
Peak values can be much larger than average, a signifi- 
cant design issue, while the low Vs during the phase when 
propulsion is absent shows that inertial effects are small 



(P < corresponds to deceleration). Each swimmer 
has its own characteristic behavior reflecting the mode 
of propulsion. In the examples shown, high power is 
needed by biflaps early in the cycle to build up speed 
that then drops steadily, tube has an even more sharply 
peaked early speed, for legs the increase in power and 
speed occurs as the flaps approach the rear, while the 
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FIG. 3: Dependence of (a) instantaneous power and (b) swim 
speed on phase angle. 



values are almost constant during the sliding portion of 
the cigar cycle. 

In conclusion, the capabilities of a variety of swim- 
mer designs have been analyzed using an MD approach 
built on particle-based swimmers and fluid. The results 
demonstrate that MD is a viable tool for modeling swim- 
ming on the microscopic scale and, more specifically, that 
swimming efficiency varies widely among swimmer body 
designs and propulsion mechanisms. Since MD is limited 
only by computational resources, it has the potential to 
play an important role in the study of microrobotic swim- 
mers, a field holding enormous industrial and therapeutic 
promise. 
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